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Abstract. Paramagnetic alignment of fractal suprathermally 
rotating grains is discussed. It is shown that if the concen- 
tration of H2 formation sites is low and resurfacing is active, 
fractal structure of grains enhances their alignment. Studying 
the influence of grain surface physics and chemistry on the 
alignment we found that there exist two critical values of grain 
sizes, and the alignment of smaller grains is expected to de- 
crease as compared to the predictions of the Purcell's theory 
(1979). One of the critical sizes is due to intensified poisoning 
of active sites, while the other is a result of a finite value of the 
imaginary part of magnetic susceptibility in the limit of high 
angular velocities. On the contrary, if active sites completely 
cover grain surface, suprathermal rotation, and therefore effi- 
cient alignment, is possible only for a limited range of grain 
sizes. 



(Spitzer & McGlynn 1979) to improve the alignment only 
marginally if high rates of resurfacing were assumed. Here we 
study physics and chemistry of resurfacing and the changes 
that a fractal structure of grains entails. The assumption of 
porous irregular grain structure, which can be described using 
a notion of fractality, corresponds to the present-day state of 
our knowledge (see Wright 1987, Whittet 1992, Williams 1993, 
Lume & Ranhola 1994). 

The structure of this paper is as follows. The concept of 
fractality in application to grains is briefly discussed in sec- 
tion 2, while the suprathermal rotation of grains due to recoils 
from H2 molecules is studied in section 3. Then in section 4, 
we address the paramagnetic mechanism of grain alignment 
and the effect of grain fractality on the measure of alignment. 
In section 5, theoretical predictions are analysed in the light of 
observational data available. 
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1. Introduction 

Although the starlight polarization by the ISM dust was 
first reported as far back as late 40s (see Hiltner 1949) it 
is still a theoretical challenge to explain the alignment of 
dust grain that causes it. The most widely accepted ex- 
planation of this phenomenon as being due to paramag- 
netic relaxation (Davis & Greenstein 1951) was later criticised 
(Jones & Spitzer 1967) as inadequate to account for observa- 
tions unless enhanced imaginary part of grain magnetic sus- 
ceptibility is assumed (see also Duley 1978, Sorrell 1994). This 
assumption became a corner stone of the Mathis (1986) the- 
ory of alignment of superparamagnetic grains. Another at- 
tempt to overcome the said drawback was made in Purcell 
(1975, 1979), where spinning of grains up to suprathermal 
velocities was considered to diminish the influence of ther- 
mal disorientation. These two modifications of the param- 
agnetic mechanism are widely believed to be the strongest 
candidates for explaining starlight polarization, though there 
are indications that Gold-type processes (Gold 1951, 1952, 
Dolginov & Mytrophanov 1976, Lazarian 1994a) are rather 
common in the ISM. 

In this paper we concentrate our attention on the Pur- 
cell's (1979) modification of the paramagnetic mechanism. 
For grains approximated by regular ellipsoids, it was shown 



2. Fractal grains 

Whatever the cause of grain formation they are expected to 
be irregular. So are the species found in the stratosphere 
(see fig. 7.3 in Williams 1993) and those formed in an arc 
(Codina-Landaberry et al. 1974). In fact, it is believed that 
ISM grains may consist of small particles jumbled together 
loosely (Mathis & Whiffen 1989). Such grains 'are the natural 
result of coagulation and disruption of grains as they cycle into 
clouds' (Mathis 1990). 

Therefore it seems natural to assume that ISM grains are 
fractal (Mandelbrot 1983). This idea was advocated in Wright 
(1987) on the basis of the long- wavelength adsorption data, and 
here we are going to study its implications for the alignment. 

To remind our reader a few facts about fractals, we use 
probably the most familiar example of fractals related to the 
length of a coastline. This length G depends on the scale of 
measurements Go as 

G ~ G T ~ d (1) 

where d is the fractal dimension, and T is the topological di- 
mension, which is 1 for the coastline. Turning from the coast- 
line to grains, we consider the surface of a well defined volume 
Vq. In this case, one may state that two grains are similar if 
for them the ratio Sy d (ro)/V^ 3 is the same, where Sf(ro) is 
the grain surface area measured with the yardstick ro . If grains 
consist of small spheres, the size of the latter determines ro. 
The analogy between a coastline and grains can be protracted 



2 



as although in both cases we do not know all the details of 
their formation we are still able to define their dimension. 

Further on for the sake of convenience, we will use I for the 
dimension of grains (a diameter for spherical species) instead of 



1/3 



An important property of fractal grains that influences 



our results is that Sf ~ l d . The actual value of d may be 
calculated within an adopted model. If the spheres of size ro 
are attached randomly one by one to a growing cluster, their 
number grows as I 2 ' 5 (Witten & Cates 1986). Assuming that 
the area of intersection is small, it is easy to see that 
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where d = 2.5 should be put for the latter case. However, if 
the spheres themselves are porous, the actual surface can be 
greater. The fractal dimension in its turn can change with the 
scale that we study. If grains are being formed from loosely 
connected PAH structures (Duley & Williams 1993), the ratio 
l/ro may reach 10 3 for grains of I ~ 3 • 10 -5 cm. Yet, we will 
use more conservative estimates in our paper. 

In short, grains being three dimensional objects can have 
surface growing faster than I 2 . In contrast, if grains are formed 
by chains of spheres sticking together, these loosely complexes 
can have dimension less than 2 (Wright 1987). Later we will 
be mostly concerned with the three dimensional clusters, as 
we believe that very loose complexes are not strong enough to 
withstand suprathermal rotation. 

3. Suprathermal rotation of fractal grains 

The concept of suprathermal rotation was first introduced in 
Purcell (1979). Various processes, e.g. variations of the ac- 
commodation and photo-emission coefficients, variations of the 
number of sites of H2 formation and the difference in cross- 
sections for left- and right-hand circularly polarized quanta can 
result in suprathermal rotation of grains. An important point 
(see Lazarian 1994b) is that the processes associated with cor- 
puscular and radiative fluxes bring about the Gold-type align- 
ment which, as a rule, has a much shorter time-scale. Therefore 
for the rest of the paper we will consider rotation induced by 
the H2 formation over the grain surface. These reactions are 
expected to proceed with a rather high efficiency: almost ev- 
ery H atom adsorbed by the surface must leave it as part of 
H2 molecule (see Williams 1988). The surface concentration of 
active sites varies for different grain species. 

We begin with considering low density distribution of ac- 
tive sites. As a rule, the sites are subdivided into two cate- 
gories: sites of physical adsorption with the energy less than 
0.5 eV and those of chemical adsorption with the energy 
of ~ 0.5 — 5 eV (Watson 1976). Although only a few sites 
of chemical binding per grain are required to account for 
the observed quantities of H2 (Hollenbach & Salpeter 1970, 
1971), some authors believe that their number is much greater 
and corresponds to the density of aH2 = 10 cm - (see 
Tielens & Allamandola 1987). 

Due to quantum tunnelling, adsorbed H atoms acquire 
high mobility (Watson 1976) and form H2 molecules both over 
external grain surface and within numerous pores of a frac- 
tal grain. The effects of these two types of reactions upon the 
grain are different. It is obvious that if an active site is lo- 
cated well within a narrow pore, the ejected H2 molecule has 
much chances to thermolise by colliding with the pore walls 



before leaving the grain. The entire energy of the reaction is 
expected to be transferred to the grain in the form of heat 
and the recoil from such H2 ejection is negligible. For active 
sites on the external grain surface, a recoil from every depart- 
ing molecule is m.H2VH2 = 2(rau2-E) 1 ^ 2 , where E is the kinetic 
energy of the ejected H2 molecule ~ 0.2 eV (Williams 1988). 
The latter value intimately depends on the chemical composi- 
tion of grain surface, and this can bring about differences in 
grain alignment for different species (Lazarian 1994c). The sit- 
uation is intermediate for grains within pores but close to the 
external surface. They are expected to experience a few colli- 
sions with pore walls before leaving the grain and the residual 
energy of the molecule depends on the details of its interaction 
with the grain surface in the course of collisions, e.g. on the ac- 
commodation coefficient. This includes first of all the portion 
of chemical energy converted into H2 molecule translational 
motion. Apart from this energy channel, there are other possi- 
bilities, e.g. the transfer of vibrational energy of a newly formed 
H2 molecule to vibrational modes of the grain. Because of the 
inherent quantum nature of this process the contributions to 
different channels depend on the spacing of the energy levels in- 
volved as well as on the nature of the interaction. The estimates 
(see Lucas & Ewing 1981) indicate that a relatively small part 
of recombination energy is transferred into translational mo- 
tion. The ejected H2 molecules can still be vibrationally 'hot', 
and this may influence their later interactions with the grain. 
Indeed, a molecule colliding with a grain loses a fraction of 
its translational energy that depends on the ratio of the De- 
bye frequency lob and the inverse time-scale of the collision loo 
(Watson 1976) 
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where Eh is the binding energy of H2 molecule, and b is 
the length scale for the repulsive interaction ~ 10 -r cm. For 
u)o <C ud, the energy transfer to the grain is ~ (E + E o )l0qL0~^ 
(Hollenbach & Salpeter 1970), which enables a molecule to col- 
lide ~ loqLO~^ times with the grain before losing its transla- 
tional energy. At the same time, one cannot exclude the possi- 
bility of converting some additional internal energy of 'hot' H2 
molecules into translational energy during these interactions 
(Watson 1976). This can make internal active sites well within 
the body of a fractal grain efficient for ejecting high speed H2 
molecules through its pores. 

For simplicity, we consider that molecules formed within 
pores up to l ac depth from the grain outer boundary re- 
tain their energy. For fractal porous grains the core of scale 
I — lac is passive in terms of contributing to the recoils from 
H2 molecules, although H2 formation is expected to proceed 
within this 'passive' core nearly with the same rate as over the 
outer surface. Among other factors, l ac depends on the pore 
geometry. For pores approximated by cylinders, it is only the 
normal component of velocity that decrease due to molecular 
collisions with the pore walls. Integrating over a Gaussian dis- 
tribution of emerging molecules with the mean energy E, it is 
possible to show that the grain momentum gained as a result 
of emission of n H2 molecules through such a pore is equal 
to n (_Emij2/37r) 1 ^ 2 and the torque is equal to n E/3. Note, 
that we considered a pore with one end closed and only the 
velocity component along the pores axis preserved. The torque 
is less if the pore is irregularly twisted. 
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It was shown in Purcell (1979) that internal dissipation of 
energy within grains, mainly due to the Barnett relaxation, 1 
suppresses rotation around any axis but the axis of the greatest 
inertia on the time-scale ~ 10 /n s, where n is the ratio of 
grain rotational energy to the equipartition energy ~ kT. Thus 
grains rotate around their major axes of inertia that we denote 
here 2-axis. Therefore the problem becomes one dimensional 
as the two other components of the torque contribute only to 
insignificant nutations. The number of H2 molecules ejected 
per second from an individual site is ~ 71 1 v\nuv~ , where 71 
is the portion of H atoms with velocity v\ and concentration 
ran adsorbed by the grain, while v~ is the number of active 
sites over the surface grain. Then the mean square of a residual 
torque is 
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which spins up grains to angular velocities limited only by fric- 
tion forces, i.e. 



Q = ([M Z ] 2 } 1/2 ^ 
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(5) 



where I z is the z component of the momentum of inertia, and 
td is the rotational damping time (Spitzer & McGlynn 1979): 



t d « 0.6i m = 0.6. 
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where t m is the time that takes a gain of density g s to collide 
with gaseous atoms of the net mass equal to that of the grain, 
ra and v\ denote the density and velocity of atoms, respectively, 
while m is the mass of an individual atom. Therefore 
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which is 

n « 1.4 • 10 
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if we assume that a grain of I = 2 • 10 - cm has v = 100 active 
sites and d = 2.5. If active sites completely cover the grain ctH2 
becomes ~ 2 • 10 15 (see Tielens & Allamandola 1987), and Q 
reduces to 1.4 • 10 s - which is comparable with the thermal 
rotational velocity of grains for T g = 90 K and T s = 10 K 



LOT « 9 • 10 



2 • 10" 



0.5(T g + iy 
50 K 



(9) 



and for grains with 1x2- 10 -6 cm, lot becomes greater than 
fi. If lac < I, thermal rotation dominates for a wide range of 
grain sizes. 2 



1 A quantitative treatment through solving Fokker-Plank equa- 
tions is given in Lazarian (1994a). 

2 The mean temperature in Eq. (9) is used because atoms ther- 
molise over grain surface after non-elastic collisions. 



Note, that fractal grains with small l ac /l ratio are likely to 
lose their suprathermal advantage. This corresponds to the sit- 
uation when the absolute majority of H2 molecules thermolise 
before leaving the grain. The H2 thermalization depends on the 
material of the grain, its size and the structure of its fractal 
maze. 

The implications of the increase of the number of active 
sites may be dramatic. It is possible to show that for a regularly 
shaped fractal grain with a uniform distribution of active sites 
chemically driven rotation corresponds to temperature 
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where T a is the gas temperature, and k is the Boltzman con- 
stant. For T a = 80 K and l ac /l = 0.3, this gives just 100 K, 
which in fact nullifies the effect of suprathermality. To preserve 
the said effect, we may allow for the change of the density of 
effective active sites per unit of grain geometric surface au2 
over patches of area eH2 The corresponding variation 
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(ctH2) 
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will be used further. We emphasise effective, since only ac- 
tive sites in pores within the depth l ac are considered. As the 
number of pores over grain surface varies, this gives a natural 
explanation for variations of ctH2- 

If we ignore the impacts from H atoms, the pressure due 
to the recoils from H2 molecules is 
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where mnnHVi is the mass flux of H atoms. Therefore the 
torque applied to the grain due to the variation of ctH2 over a 
patch CH2 is proportional to 

1 2 7i 
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where it is taken into account that the number of ch2 patches 
is 

(^) • The 

mean angular velocity of grains is 

which gives ~ 5 • 10 5 s -1 for ran = ra, l ac = I and ( c h2~®h2 ^ _ 
1 

300 ' 

Assuming that &H2 is proportional to l 1 ^ 2 , it is easy to see 
that Q ~ I -3 / 2 and therefore at sizes ~ 5 • 10 -6 cm the ther- 
mal effects become dominant. Similarly, for large I, a greater 



3 Naturally, a more complete description of surface variations 
may be given by employing a distribution of scale lengths, i.e. 
a spectrum of spatial frequencies. In our model, we avoid these 
complications. 
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portion of molecules thermolise before leaving deep pores of 
grains and rotation cease to be suprathermal. 

Therefore if fractal grains are completely covered by active 
sites, they can be suprathermal due to H2 formation within 
a limited range of sizes [ii,^]; for grains larger than h and 
smaller than h thermal effects dominate. Whether or not these 
ranges overlap as well as the precise values of h and h de- 
pend on the grain composition and structure. We intend to 
address this problem in our future research. Note that it is due 
to grain fractality that the upper limit h appears; for solid 
smooth grains covered with active sites only the lower limit 
l\ ~ 10 - cm is expected. 

4. Paramagnetic mechanism 
4-1. Factors of alignment 

A paramagnetic grain rotating in the external magnetic field B 
is subjected to dissipative torques (Davis & Greenstein 1951 , 
see also Jones & Spitzer 1967, 

Purcell & Spitzer 1971, Roberge et al. 1993). In the grain co- 
ordinate system, the component of B perpendicular to the z 
axis appears as a rotating field. If the angle between B and z 
is 9, the magnitude of the said component is Bsm9. The mag- 
netic moment M produced by the imaginary part of the mag- 
netic susceptibility x" i s perpendicular to both z and B and 
has the magnitude Vx" B sind . The retarding torque M x B 
causes rotation of the vector of angular momentum J = Iii 
towards the direction of B, namely, 



d0_ _ V X "B 2 sing 

dt ~ TTo 



(16) 



where it is taken into account that the rotation is around the 
axis of maximal inertia (i.e. / = I z ). 

For low frequencies, x" i s given by (Spitzer 1978) 
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which holds up to frequencies ui c = 4 x 10 11 N e /i 2 /k, where 
fi = 0.9 xlO 20 erg G -1 is the electron magnetic moment and N e 
is the number of unpaired electron spins per cm . For a grains 
with one paramagnetic ion per hundred atoms, lo c ~ 10 s -1 . 
If 3> Lo c , which may be a case for suprathermally rotating 
grains, x" becomes independent of Q and approaches the static 
susceptibility. Therefore at frequencies 10 — 10 s - , x" is ex- 
pected to come to saturation and not to depend on Q. Note, 
that the dependence Tg is preserved for these saturation val- 
ues. 

In short, large angular velocities of grains increase the char- 
acteristic time of paramagnetic relaxation as x" comes to its 
asymptotic values (compare Whittet 1992). The characteristic 
time of this relaxation is 
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where Ts is the temperature of grain material, q is its den- 
sity, q is the factor of imaginary part of magnetic susceptibility 
enhancement (Jones & Spitzer 1967), and ii'(fi) is a function 
with the following asymptotics 
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For our model, we will consider that transition from one asymp- 
totic to another is at Q = co c . 

For the original Davis-Greenstein process, t r should be 
compared with the time that takes a grain of mass m to collide 
with gas atoms of the same total mass 



2 • 10 1 



I 



2 • 10 5 cm 



g cm" 



1.3 • 10 5 cm s" 1 



20 cm" 



10" 



(20) 



where the velocity v a corresponds to that of H atoms at T fa 
80 K. Note, that in many cases the velocity v a is dominated by 
Alfvenic phenomena (see Aron & Max 1975, Elmegreen 1992, 
Myers 1987, 1994, Lazarian 1992, 1993, 1994a), the fact that 
makes t m even smaller. The ratio of t m to t r 
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determines the mean value of precession angle fi and therefore 
the measure of grain alignment is 
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According to J. Mathis (private communication), the low value 
of B in Eq. (21) should be treated as standard and 'there might 
be as much as a fractor 10 or more increase in <5i for the B 
term (see Boulares & Cox 1990, fig. 3b)'. It is possible to show 
that for the Davis-Greenstein process a takes the following 
form (Aannestad & Greenberg 1983) 
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The corresponding plot in Fig. 1 indicates that it is only 
through an appreciable increase of q, as compared to the 
usual materials, (see Jones & Spitzer 1967) that one can rec- 
oncile the mechanism with observations. Although not im- 
possible (see Sorrell 1994), the above conjecture needs addi- 
tional assumptions (Mathis 1986) to explain the polarization- 
wavelength dependence (see Wilking et al. 1982). There are 
also indications that the modified theory does not fit the UV 
observations (see Clayton et al. 1992). 

Thus, the suggestion that a better alignment can be de- 
livered by suprathermally rotating grains (Purcell 1975, 1979) 
was met enthusiastically in the community (see Martin 1978, 
Johnson 1982). Indeed, if grain rotational energy is greater 
than k times any temperature in the system, chaotic bombard- 
ment of ambient gas cannot alter the regular rotation. Thus, if 
the time associated with a regular torque (see section 3) is ti, 
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Fig. 1. The measure of the Davis-Greenstein alignment as a function 
of grains size measured in 10 — 5 cm units for different factors q. 



it is this time that should be compared with t r . In general, II 
may be much greater than t m and therefore it gives a chance 
of improving alignment without forcing t r down. It was argued 
in Purcell (1979) that the quasi-empirical formula 



i 0.6 




t x = 1.3(ii, +0.6i m ) 



(25) 



is a good approximation within the accuracy of 20% for the 
'mean time back to crossover' and this should be compared 
with t r . Indeed, if II <C t m , the applied torques are effectively 
chaotic and t x should be ~ t m - In the opposite case (ti ^ 
t m ), which is much more interesting in the view of the above 
discussion, t x should be of the order of ti which is also satisfied. 
Thus, although Eq. (25) is not a precise theoretical result and 
was only tested for 0.1< ti/tm < 10 (Purcell 1979) we will use 
it for the whole range of ti that we encounter further in this 
paper. 

In was shown in Spitzer & McGlynn (1979) that if the 
disorientation is not complete during a crossover the following 
parameter 8 e fr enters the theory: 



F>1 
Jx- F < 1 

Ft r ' 
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where F is the 'disorientation parameter', which is discussed 
further. The mean alignment measure [see Eq. (22)] may be 
expressed using 8 e g (Aannestad & Greenberg 1983) 



(27) 



a(8) exp < — - — > d<5 

b eff 



where cr(8) is given by (Purcell 1979) 

3 1- (exp 8 - 1)" 1/2 arctan(exp 8 - 1) 1/2 1 
a( ' ~ 2 [1 - exp(-<5)] 2 

This measure is shown in Fig. 2 as a function of 8 e g . 
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4.2. Spin-up for grains 

It follows from the above that to account for the alignment 
when t r is given by Eq. (18), one needs to have ti 3> t m - The 
difficulty on this way is the finite life time of active sites over 
grain surface. 



Fig. 2. The measure of suprathermal alignment as a function of 6 e 



Active sites of H2 formation are believed to be associated 
with places of chemical binding. H atoms physically adsorbed 
by a grain (the binding energy < 0.1 eV) migrate over its 
surface until they meet either a vacant site for chemisorption 
or a trapped chemisorbed H atom and form a molecule (see 
Tielens & Allamandola 1987). 

A process that can limit the life of active sites is 
resurfacing due to accretion of a new monomolecular layer 
(Spitzer & McGlynn 1979). The accretion rate is proportional 
to the incoming atomic flux times overall grain cross-section 
S(l) ~ 7rl . Therefore a grain of mass m(l) grows as 



dm(Q 
At 
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S(l)v a n a m a ^ a 



(29) 



where f a is the adsorption probability, a coefficient with the 
range of values < f a < 1 for atoms of density n a , velocity v a 
and mass m a - Assuming that accreting atoms uniformly cover 
the entire physical surface of grains and substituting dm(l) = 
Q,Sf(l)dl, where Sf(l) is given by Eq. (2), one obtains for the 
accreting time: 
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where Ah is the thickness of an adsorbed layer. If vh and 
v a are due to thermal motions, then = (7^") ^ an d 

t a becomes proportional to ^ • However, if non-thermal 

Alfvenic motions dominate, vh becomes equal to v a and t a is 
less or proportional to (^f L ) ■ The ratio (l/ro ) d ~ 2 , where ro is 
the size of an elementary sphere our model grain is composed 
of, appears due to assumed high mobility of accreted species. 
Whereas if poisoning atoms hydrogenate on their arrival, the 
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products like H2O are essentially immobile over grain surface. 
In such a case, accretion is active only over external grain sur- 
face. However, the life time of active sites within grain pores is 
expected to increase as compared with that discussed in Spitzer 
& McGlynn (1979). 

We believe that oxygen is the most important agent poi- 
soning the active sites. Other chemical elements either form 
compositions that desorb more easily than H2O or have a 
negligible abundance in the gaseous phase. As all heavy el- 
ements oxygen is usually depleted in the gaseous phase. For 
instance, its depletion towards £ Persei is -0.5, which gives 



no/ni 



10" 



(Cardelli et al. 1991). Assuming that a layer 



of A/t PS 3.7 • 10 8 cm effectively covers active sites over grain 
surface (Aannestad & Greenberg 1983), one obtains 
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1.85 • 10" 
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Therefore for thermal accretion one has 
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considering the accretion of oxygen atoms. The estimate for 
t a can be further increased if one assumes that f a < 1 and 
the desorption is present. The value of f a ~ 0.1 was used in 
Spitzer (1978, p. 208) for illustrative purposes while values 0.02 
and 0.2 were reported for sticking probabilities of C over sil- 
icate and graphite covered with a monolayer of H2O surfaces 
(see Williams 1993). It was shown in Jones & Williams (1984) 
that the effective sticking probability varies as mantles begin 
to be deposited and as a consequence the time of establish- 
ing monolayer coverage may be increased by a factor of about 
three. 

The products of grain surface chemistry are also subjected 
to desorption which increases the time of resurfacing. Accord- 
ing to Jones & Williams (1984) the probability of retaining the 
products OH and H2O on the grains surface is not less than 
0.7 in dark quiescent clouds, whereas NH may be predomi- 
nantly released on formation (see Wagenblast et al. 1993). Im- 
mediately after formation the reaction products are subjected 
to photodesorption. This process can be driven not only by 
UV quanta (see Johnson 1982), but also by water ice 3 /tin 
absorption (Williams et al. 1992). Note, that 3 /tin radiation 
penetrates much further into dark clouds. The energy release 
associated with H2 formation may be another cause of des- 
orbing weakly bound molecules such as CO in the vicinity of 
active sites (Duley & Williams 1993). This means that active 
sites are able to control to some extent resurfacing in their 
vicinity. 

In brief, we may expect to obtain t a of the order of 
10 tm(l/ro) d ~ 2 using rather conservative estimates for the 
above parameters. 

The assumption of a uniform coverage of grain surface by 
heavy atoms may not be true. Indeed, the motion of heavy 
atoms on the surface is inhibited by the irregularities of the 
surface potential that result from numerous surface defects and 
the discrete position of the lattice molecules. The exact values 
for the corresponding potential minima are not known for the 
interstellar grains and we have to relay on the laboratory data 
for the materials which probably have many fewer defects than 
those in the ISM. Atoms trapped in lattice defects may enter 
a potential minimum with an energy comparable to that of 



chemisorption and lose their mobility for the conditions we dis- 
cuss here. For physically adsorbed atoms to move between ad- 
jacent minima (a few Angstrom) a barrier of Eb height should 
be surmounted. Eb is typically around one third of the total 
binding energy (although see Jaycock & Parfitt 1986) which is 
greater than 800 K for O, N and C atoms (Watson 1976). The 
probability of quantum tunnelling for heavy atoms is negli- 
gible and therefore their mobility is due to thermal hopping. 
The characteristic time for a particle to hop over a barrier is 
(Watson 1976) 



th 



1 (E B 
— exp 1 



\kTs 



(33) 



where vo ~ 10 12 s -1 is a characteristic vibration frequency for 
the adsorbed particle. For T s = 7 K, this gives th ~ 10 5 s. 
Since this is a random walk process, a time 



^ M2exp Gt) 



(34) 



is required for a particle to get to a point M lattice spacing 
away. This makes the diffusion time for a heavy atom inside 
an irregular fractal structure of a grain of 10 -5 cm compara- 
ble with the times of grain-grain collisions, i.e. greater than t r . 
Such situation is expected for a low abundance of H atoms or 
when reaction between H atoms and 'poisoners' are strongly 
inhibited. On the contrary, as long as H atoms have a frac- 
tional abundance greater than about 10 -3 , there are enough H 
atoms to hydrogenate, essentially instantaneously, every heavy 
atom that is accreted (see Millar 1993). The products of the 
reactions are expected either to be ejected as OH, CH, CH2, 
CH3, NH, N H2 or remain on the grain to produce fully satu- 
rated forms like H2O, CHu, NH3. These molecules are highly 
immobile and their diffusion inside pores is negligible. There- 
fore the probability of covering active sites even relatively close 
to the outer surface of a grain may be substantially reduced as 
compared to resurfacing of the outer parts of the grain. How- 
ever, to dominate grain dynamics the the number of active sites 
within the grain up to the depth l ac should exceed the number 
of active sites over the outer surface. 

Our discussion reveals that for fractal grains simple esti- 
mates for tL given in Spitzer & McGlynn (1979) are no more 
applicable. In fact, tL apart from the subtle details of grain 
surface chemistry may depend on grain fractal dimension as 
well as on I and l ac . To illustrate this phenomenon, we may 
assume 



tL = zut„ 



I 



(35) 



where vj is a parameter of the model, which is chosen 10 in 
accordance with the considerations above. Note, that our con- 
jecture for ti is not applicable even qualitatively for small I. 
Indeed, if the time of grain collision with H atom 



tc 



1 



■ktihvhI 2 



3 • W 



2 • 10" 



I 



20 cm" 



80 K 



(36) 



becomes less than the average time for an oxygen 
atom to scan the grain surface with As sites (see 
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Tielens & Allamandola 1987): 



Asvo 1 ( E B 

-^ exp Km 



3 • 10° 



x exp 



2 • 10" 5 cm) lr 
f E B /k\ fWK 



\ 800 K 



(37) 



oxygen atoms poison sites of H2 formation before a 
monomolecular layer is formed. For grain surfaces cov- 
ered with a few monolayers of ice (see Zhang & Buch 1990, 
Buch & Zhang 1991) this may mean that oxygen occupies 
places which correspond to the potential minima, e.g. places of 
most probable H2 formation. The critical size is 



10" 



x exp 



d-2 / 

— ( 20 cm"' 

V ~ 

E B /k\ [10 K 



4 ( T * 
V80 K 



800 K ) V T 3 



(38) 



which gives 0.6 • 10 -5 cm for l cr /ro = 100 and d = 2.5. Nat- 
urally, grains with a few active sites are much more l cr sen- 
sitive than those with high density of active sites. Indeed, 
if a grain with I = 0.5 • 10 - cm has only 10 active sites 
(Hollenbach & Salpeter 1971, Purcell 1979), it needs to collide 
just with a hundred O atoms to cover them. This number be- 
comes more than 3-10 if the grain is completely covered by 
the sites. For I less than l cr we expect a drastic decrease in 
II, which is inversely proportional to the fraction of surface 
covered by active sites of H2 formation and th becomes much 
less than t m . As soon as ij, < < m , one has to use th instead 
of td in Eq. (6). Thus the 'regular' Q falls much lower than 
'thermal' lot [see Eq. (9)] consequently the rotation ceize to 
be suprathermal. The original Davis-Greenstein theory should 
be applied to obtain the measure of alignment [see Eqs (23) 
and (14)]. 

To summarise, it has been shown that th depends on subtle 
processes of grain surface chemistry and for grains greater than 
lcr, th is expected to be greater than the estimate in Spitzer 
& McGlynn (1979) if grains are fractal. 

4-3. 'Crossovers' 

As the z component of the torque acting on each suprather- 
mally rotating grain (see section 3) changes its sign, the grain 
will be braked and then will start rotating around the axis 
which, in general, is not parallel to the initail axis of rota- 
tion. However, it was shown in Spitzer & McGlynn (1979) that 
a partial preservation of orientation during such crossovers is 
possible. This provides for the disorientation factor (see Ap- 
pendix) 



F«0, 



rriH2VH2 



g cm s ± ) \ 3 g cm 
I \ 3 ( n H 



v H 
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1 \ — 3 — / / 

l ac \ I l ac 
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I 



-1/6 



(39) 



where l ac is assumed greater than To. The factor (an2/4 • 
10 10 cm" 2 ) 1/3 7rl 2 (l ac /r ) (d " 2)/3 is the number of active sites 
contributing to the spin-up. The disorientation grows with this 
number and decreases with the fraction of surface involved in 
'torque generating', i.e. (^p) . 

4-4- Grains completely covered by active sites 

If a suprathermally rotating grain is completely covered by 
active sites and the uncompensated torque is due to the grain 
shape, this shape has to be substantially altered, e.g. due to 
grain-grain collisions, to change the direction of the torque. 
The corresponding time is, as a rule, greater than t r . 

At the same time, we have seen in section 3 that the degree 
of suprathermality is low for such grains. Therefore fluctuations 
due to the tail of the Maxwellian distribution are likely to in- 
fluence them. This influence is expected to increase for large 
grains as they lose suprathermality due to (l a c/l) coefficient 
and for small grains due to a rapid growth of lot with the 
decrease of grain size. For some materials and fractal dimen- 
sions these sizes can overlap, and such grains cannot exhibit 
suprathermal behaviour. Their alignment should be described 
by the Davis-Greenstein theory [see Eqs (23) and (10)]. 

If the spatial variations of H2 formation are due to physico- 
chemical properties of grain surface , crossovers are expected 
as a result of resurfacing and our earlier estimates for ti, are 
applicable, provided that grains are suprathermal. 

4.5. Grain with a low density of active sites 

If there are only a few active sites over grain surface, the ac- 
cretion of a monomolecular layer alters substantially grain ro- 
tation. For F < 1 and I > l cr , one obtains 
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(40) 



If ti, is given by Eq. (35), it is possible to plot the alingment 
measure as a function of I (see Fig. 3). 

It is worth noting that the fractal dimension can vary with 
grain composition, size and physical conditions around. For 
instance, grain pores can be covered by ice as a result of mantle 
growth. If ti, ^> t m , the alignment measure for large grains 
decreases less rapidly because of partial preservation of the 
angular momentum direction during crossovers which results in 
F < 1 (see Fig. 3). The increase of the said measure of claimed 
in Spitzer & McGlynn (1979) is possible for th < t m . The 
discontinuity of the derivative that is obvious from Fig. 3 is an 
artifact of the simplified model adopted. Indeed, even for F > 1 
there should be a residual preservation of alignment, while in 
our model a complete disorientation after the threshold F = 1 
is assumed. In fact, one should bear in mind that Fig. 3 is 
obtained within a model which adopts a conjecture for th given 
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by Eq. (35). Further research should be directed to obtaining a 
more realistic analytical approximation for the time of spin-up 



Our calculations in section 3 have shown that if the density 
of active sites over grain surface is high, such grains are likely 
to lose the suprathermal advantage if their sizes are less than 
2 • 10 -6 cm. At the same time, our calculations in section 4 
indicate that if there are only a few sites of H2 formation per 
grain, poisoning mainly by O atoms intensifies for I less than 



l cr 



6 • 10 cm. Therefore whatever the density of active 



°- 4 

a 



d=2.8 
d=2.5 
d=2.2 



sites, grains less than a particular size are likely to be aligned 
only marginally. 

The critical size l cr can increase due to the inverse depen- 
dence of Ts on I (Greenberg 1971): 



T s (l) « T s 



2 • 10" 



I 



K 



(41) 



1 



where Ts is the temperature of a grain of size 2 • 10 - cm. This 
is a result of suppressed thermal emission from small grains in 
infrared. Although the power law index is just 0.2, this depen- 
dence may be essential in the view of exponential dependence 
of lcr on temperature [see Eq. (38)]. 

Another important factor is the difference between the 
imaginary parts of magnetic susceptibilities of slowly and 
rapidly rotating grains. Since the relaxation time for Q > 
10 9 s -1 is proportional to Q (see section 3), 8\ ~ (Ifi) -1 . As Q 
increases as | - ( 2 + d ) with the decrease of I, Si should decrease at 
least as I -1 . Thus, we expect to see a decrease of alignment for 
I smaller than a critical value. For a grain with a few hundreds 
of active sites, this corresponds to I < 10 -5 cm. 

Note, that even for thermally driven rotation given by 
Eq. (9) we expect the imaginary part of magnetic suscepti- 
Tig. 3. The alignment measure as a function of grain size I in units bility to come to the saturation for I ~ 10" 6 cm. Therefore 
of 2-10" 5 cmforro = 2 • 10" 6 cm and different fractal dimensions ^ giyen by Eq ( 2 1) decreases for smaller I as I 3 ' 2 . In fact, 

d. We intentionally start from / = 10 cm, as for smaller grains , n n , r r ■ j 

suprathermal rotation tor some range ot gram sizes and some 

we expect the decrease 01 the alignment measure, which depends on . . i . i . r . A 

, . , , . , materials can suppress alignment instead ot increasing it. As 

gram physics and chemistry. . . . — 1 r \n 

the rate of accretion is proportional to I [see Eq. (30)], it is 

easy to show using Eqs (8) and (18) that 8 e g given by Eq. (26) 
decreases as l (d ~ 2 ^ 2 for Q > LOcr- The combination of this effect 
with enhanced poisoning for small grains makes the decrease of 
the alignment measure with I much more rapid. The variations 
in concentration of paramagnetic ions, otherwise marginally 
important, 5 may result in substantial changes in the sizes of 
aligned suprathermally rotating grains. Similarly, the decrease 
of atomic hydrogen concentration which mitigates the degree 
of suprathermality can in some cases increase alignment for a 
particular range of grain sizes. 

Other factors that can contribute to a lower degree of po- 
larization from small grains may be their shape and the de- 
creased efficiency of H2 formation over their surfaces. As for 
the shape, it is impossible to exclude that small grains become 
more spherical and the alignment of their angular momenta 
does not affect the overall polarization. As for the formation 
of H2, it is inefficient on PAH surfaces (see Tielens 1993). This 
property may be shared by some carbonaceous grains if, ac- 
cording to Duley & Williams (1988), they consist of loosely 
connected PAH structures. 



5. Theory & observations 

5.1. Small and large grains 

Using the MNR grain model (Mathis et al.1977), it is possi- 
ble to show (Mathis 1979) that a good fit of the polarization 
curve may be obtained if grains with radii above a certain 
threshold are aligned [although see Johnson (1982), where the 
observational evidence for better alignment of small grains is 
claimed]. The above dependence was accounted for in Mathis 
(1986) by assuming that only those grains that contain at 
least one superferromagnetic cluster are being aligned. Within 
the adopted model, the probability of having such a cluster is 
greater for larger grains. Alternatively, for the Gold-type align- 
ment (Lazarian 1994a), it is possible to show that only grains 
larger than 10 -5 cm are sufficiently inertial to be aligned un- 
der Alfvenic perturbations in the typical ISM conditions. Here 
we are interested whether alignment of suprathermally rotat- 
ing fractal grains can reproduce this type of behaviour for the 
polarization curve, and thus need to discuss the alignment of 
large and small grains separately. 



4 Grains are expected to be aligned with their long axis per- 
pendicular to magnetic field lines, which corresponds to obser- 
vations (Lazarian 1994a). 



It is stated in Jones & Spitzer (1967) that in the limit of low 
frequency, the imaginary part of the magnetic sucseptibility is 
independent of the magnitude of individual magnetic moments. 
Therefore in this limit, nuclear paramagnetism should be as ef- 
fective as electronic paramagnetism in relaxing the interstellar 
grains. 
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To summarise, for suprathermally rotating grains the de- 
gree of alignment is expected to be low for small grains. 

A quantitative comparison of theory and observations for 
polarisation produced by dust with a power-law size distribu- 
tion and with the grain elongation varying with size will be 
given elsewhere. 

5.2. Silicate and carbonaceous grains 

Excess in polarization associated with 9.7 /tin silicate feature 
demonstrates that silicate grains are aligned in both low and 
high density regions of the ISM (see Whittet 1992). The gen- 
eral form of the polarization curve with its smooth peak in 
visible and steady decline in infrared is difficult to reproduce 
using metals or other strong absorbers such as graphite. How- 
ever, it seems that carbonaceous grains discussed in Duley 
(1993) can be aligned. These grains have a refractive index 
m = n + ik > 1 and have a porous non-absorbing matrix 
with a dilute dispersion of absorbing molecules. It was shown 
in Duley (1993) that a polyethylene matrix should not con- 
tain more than 4 • 10 19 cm -3 molecules of coronene C24-H24 
for k be less than 0.1. Such a solid satisfies the requirement 
m ~ 1 — 1.7 and allows efficient utilisation of carbon to pro- 
duce extinction. At the same time, fluffy carbonaceous struc- 
tures which may be partially graphitized by exposure to the 
interstellar radiation field are unlikely to be aligned by the 
Purcell's mechanism (1975, 1979) as the ejected H2 molecules 
are expected to have little kinetic energy. We argue in our next 
paper (Lazarian 1994c that due to peculiarities of H2 forma- 
tion over grain surface (Duley & Williams 1986, 1993) grains 
with surfaces of H2O ice, defected silicate or polymeric car- 
bonaceous material are likely to exibit an enhanced alignment 
due to suprathermal rotation, while those of aromatic carbona- 
ceous material or graphite are not. 

5.3. Cold and warm grains 

We have already seen that even small changes of temperature 
can substantially alter grain alignment. In fact, an increase of 
temperature over a few degrees considerably increases l cr due 
its exponential dependence on temperature. As a result, poi- 
soning atoms become sufficiently mobile to attack active sites 
directly and therefore the time of changing position of active 
sites becomes less than the time of accreting a monomolecular 
layer. 

If the temperature is further increased, the H2 formation 
becomes suppressed. For a reaction between physically and 
chemically adsorbed H atoms to occur, an activation barrier 
should be overcome. A critical temperature T cr can be defined 
for which the probability of reaction with a chemisorbed atom 
is equal to the rate of evaporation: 

Tcr = , . f" . (42) 
k m(Tyi'o) 

where r r is the reaction time-scale ~ 10 -6 s. Estimates pro- 
vide T cr ~ 50 K for H2 formation over a silicate surface 
(Tielens & Allamandola 1987). This means that, for instance, 
the aligned warm dust in the core regions of OMC-1 (see 
Chrysostomou et al. 1994) can not rotate suprathermally due 
to H2 formation. Therefore it is unlikely that the suprathermal 
paramagnetic alignment of warm (Ts > 50 K) dust be efficient. 
Thus it either the Gold-type alignment (see Lazarian 1994a) 



or Davis-Greenstein mechanism that drives the alignment of 
warm dust. 

5.4- Photodesorption 

Throughout our discussion we assumed high values of grain 
resurfacing and therefore high rates of changing locations 
of active sites. This may not always be the case with 
grains in diffuse clouds. For instance, some observations (see 
Tanaka et al. 1990) testify that ice mantles occur only to- 
wards sources with visual extinction A v > 10 (although 
see Whittet et al. 1988). This may correspond to cleaning 
grain surfaces by photodesorption, and therefore the existence 
time for Purcell 'rockets' may be greater than that given by 
Eq. (30). Indeed, the strongest poisoners for active sites are 
oxygen, nitrogen and carbon as they are usually adsorbed into 
monomolecular layer in time much shorter than t r . However, 
H atoms are expected to tunnel beyond the activation barrier 
for energetically favourable reactions of H2O, NH3 and CHu 
formation. These molecules are likely to be cleaned from grain 
surfaces by photodesorption which has a characteristic time 
~ 10 s for an unshielded regions of diffuse clouds. This time 
is usually less than t r , i.e. a better alignment may be achieved 
in comparison with the predictions we have obtained in the 
present paper. 

However, we do not have conclusive observational evidence 
in favour of the picture discussed above. On the contrary, elab- 
orate modelling of dust alignment on the line of sight to the 
Becklin-Neugebauer (BN) object in Lee & Draine (1985) seems 
to testify that the core mantle grains tend to be well aligned. 
An evidence that the degree of polarization is correlated with 
the strength of 3.1 /tin water ice band follows from observa- 
tions (see Joyce & Simon 1982) although this is by no means 
a universal tendency (Dyck & Lonsdale 1981). This seems to 
contradict the notion that keeping surfaces of grains clean con- 
siderably improves their alignment (compare Johnson 1982). 
At the same time, we probably do not observe significant deple- 
tion of elements heavier than He in clouds (see Whittet 1992) 
to attribute the decrease of poisoning to the fact that all the 
'troublesome poisoners' are adsorbed by grains. Therefore the 
possibility of the prolonged spin-up due to the fractal structure 
of grains may be valuable for explaining observations. 

At the same time, indications that alignment is system- 
atically less efficient towards stars with higher extinction 
(Tamura et al. 1987, Whittet et al. 1994) as well as low degree 
of polarization in far infrared reported in Goodman (1994) may 
testify that the alignment decreases as major part of atomic 
hydrogen is converted into molecular form. 

6. Conclusions 

The main conclusion of this study is that the fractal structure 
of grains influences the efficiency of paramagnetic alignment. 
The rotation of grains with a developed fractal surface be- 
comes less vulnerable to the resurfacing induced by accretion 
and this improves the alignment. At the same time, fractality 
can inhibit alignment if the surface coverage by active sites 
is high and the majority of H2 molecules thermolise within 
fractal pores before leaving grains. Our study reveals a non- 
monotonous dependence of the alignment measure on grain 
size. The details of the said dependence are influenced by bind- 
ing energy, the number of active sites, grain composition and 
their temperature. For a low density distribution of active sites 
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there exist a critical size caused by active poisoning and, for 
grains less than this size, the alignment is inhibited due to 
frequent crossovers. Yet another 'low' critical size is shown to 
exist due to a limited value of the imaginary part of mag- 
netic susceptibility. This weakens the alignment of rapidly ro- 
tating small grains. Considering a high density of active sites, 
we found that grains both less and greater than particular sizes 
loose their suprathermality, and the original Davis-Greenstein 
theory should be applied to describe their alignment. In both 
cases, the alignment is temperature dependent and for warm 
grains (Ts > 50 K) the 'chemically driven' suprathermal rota- 
tion is unlikely. A further study of alignment for different ISM 
regions should give an insight into complex processes of grain 
physics and chemistry. 
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A. Disorientation parameter for fractal grains 

A partial preservation of alignment during crossovers can be 
described using parameter F introduced in Spitzer & McGlynn 
(1979) as 



cos(xt) = exp(-_F) 



(Al) 



where \T is the angle of deviation of the angular momentum. 
F can also be expressed as an integral over the mean squared 
deviations due to elementary torques 



F=\\ N((A Xi y)dt 



(A2) 



where N is the frequency of the torque events. Calculations in 
Spitzer & McGlynn (1979) give 



tt | ((Ai z ) 2 ) / 3((AJ ± ) 2 ) 
+ 2((AJ Z ) 2 } 



4 1 |Jl|(AJ z 



(A3) 



in the approximation of the constant angular momentum com- 
ponent that is perpendicular to the axis of the greatest inertia. 
Note, that A J z and A Jj_ are elementary changes of the angu- 
lar momentum in an individual torque event, and the angular 
brackets denote averaging over distributions of A J z and Jj_. 

The squared dispersion of Jj_ at t = is given in Spitzer 
& McGlynn (1979, eq. 42): 

<Ji}(0) « (^j 1/3 T(l)(A a N) 1 ^(AJ z )- 2/3 I 2 J 3 

x <(AJ ± ) 2 } (A4) 

where T(x) is the Gamma function, and the coefficient of the 
Barnett relaxation (Purcell 1979, eq. 46) 



A a = 6 x 10 17 ( 



s 



.2 x 10" 5 cmj ^3 g cm" 3 J 15 K 



-l 



(A5) 



is assumed. 

The frequency of the torque events is proportional to 
the frequency of grain collisions with H atoms multiplied by 



(Lc/l) d , i.e. 



N 



5 x 10" 



1>ac 



7i 
0.2 



10- 5 cm 



20 cm- 



l 



10 5 cm s 1 



(A6) 



To calculate F from Eq. (A3), one needs to know the average 
of J^(0). Note, that (Ji(0)) is a sum o (J 2 (0)} and (Jy(0)) 
and due to the symmetry inherent to the problem 



« 2 H2 = (J 2 J0)) = (J 2 (0)} = 0.5(Ji(0)} 
For a Gaussian distribution, one gets 

1 _ 1 f f + ca 1 

(Ji(0)} " 27ra ff2 

x exp 



(A7) 



' x 2 + y 2 



2 i 2 

^ + y 



dxdy 



which can be calculated in polar coordinates to give 

1 pK 1 t/tT 

(Jl(O)} = V 2~ 2 ~ 



(Ji(o)y' 2 



(A8) 



(A9) 



Similarly the Gaussian averages can be found for [AJ Z ] -2 ^ 3 , 
where the square brackets [. . .] denote averaging over the sur- 
face of an individual grain. The dispersion of [A J z ] which we 
will denote (72 is inversely proportional to the square root of the 
number of active sites v ac over the grain surface and depends 
on grain geometry and distribution of active sites in respect to 
the 2-axis. Therefore 



-2/3\ 



1 



(T 2 V27T 



K 2/3 



exp 



1 x/ 
~2^l 



dx 



(A10) 



which gives (see Gradshtein & Ryzhik 1965, [3.462(1)] and 
[9.241(1)]) 



-2/3\ 




(All) 



where D p (x) is the parabolic cylinder function and T(x) is the 
Gamma function. 
Thus 



F 



((AJ Z 



(A a 7V) 1 / 6 /i /3 ((AJ z ) 2 ) 1 / 3 ((AJ ± ) 2 } 1 / 2 
i I 3((AJx) 2 ) A 

2((aj z ) 2 ); 

where the numerical factor is equal to 



3 x -l/6_ 



-1/2 



(4/3) 



e 2 x 



-2/3 



dx 



Assuming ((AJ Z 



((AJ ± 



one obtains 



{(A.h) 2 ) 1 ' 2 



F = h- 

2 (A a 7V) 1 / 6 /i /3 ((AJ z ) 2 } 1 / 3 



(A12) 



(A13) 



(A14) 
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Substituting ((AJ Z ) 2 } = ^■m 2 H2 v 2 H2 as well as expressions for 
A a , N and ((AJ Z ) 2 } into Eq. (A14), one gets 

\ 1 / 3 / N 1/3 



-1/2 



50 



I \~ 2 ( n H \~ 1 f 6 



3gcm- 3 y \2 x 10" 5 cm/ V20 cm" 3 

VH \- 1/6 (ll \- 1/6 flac\- d/6 



1.3 x 10 5 cm s- 1 / V0.27 V I 



(A15) 



The number of active sites v ac over the grain surface up to the 
depth lac may be assumed to be proportional to the physical 
surface of the grain S g (l). Therefore for a fractal grain one 
obtains 



0.8 



1/2 / l \ " 4 / 3 



3 gcm~ 3 / \2 x 10- 5 cm 



/ 71 \ / l ac \ / *a 



20 cm" 3 / V0.27 VI/ V r 



1.3 x 10 5 cm s- 1 / V4 x 10 10 



(A16) 



which shows that F depends on the grain fractal dimension. 
If lac is much greater than ro and d > 2, F increases which 
means that a partial preservation of alignment during crossover 
becomes marginal for small grains. However, as the number of 
active sites over the depth l ac decreases, the preservation of 
alignment becomes more important. 
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